Noise models for superoperators in the chord representation 
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We study many-qubit generalizations of quantum noise channels that can be written as an incoherent sum of 
translations in phase space, for which the chord representation results specially useful. Physical description in 
terms of the spectral properties of the superoperator and the action in phase space are provided. A very natural 
description of decoherence leading to a preferred basis is achieved with diffusion along a phase space line. The 
numerical advantages of using the chord representation are illustrated in the case of coarse-graining noise. 
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I. INTRODUCTION 

The precise manipulation of coherent quantum processes 
is the ultimate goal at the basis of quantum information the- 
ory. However experimental quantum systems cannot be com- 
pletely isolated and therefore unwanted interactions pose strict 
limits on its practical factibility. 

Therefore understanding the effects of certain kinds of 
noise on algorithms and how to correct the errors they pro- 
duce is a major subject of interest. 

The purpose of this paper is to show how some simple noise 
models of interest in quantum information theory can be for- 
mulated in a phase space language. In doing so we reinterpret 
the depolarizing and the phase damping channels on many 
qubits as a special kind of diffusion in phase space, study their 
spectral properties and display their effects on selected pure 
initial states. 

The Weyl representation, and its associated Wigner func- 
tion, have been successfully used to understand many aspects 
of quantum mechanics 1 1], especially its classical limit |y,|3|]. 
In its discrete version, it has also been used to study the prop- 
erties of quantum maps and, recently, to analyze many-qubit 
quantum algorithms |4|. In phase space the Wigner distri- 
bution of a state displays both its classical properties, in the 
smooth part of the distribution, as well as its quantum prop- 
erties, in the form of highly coiTelated oscillatory structures 
with sub-Planck scale sizes |5]. To reveal these structures 
ocurring at widely different scales it is then natural to con- 
sider the Fourier transform of the Wigner function. This pro- 
cedure leads to a different, and alternative, representation of 
quantum mechanics in phase space. This representation, when 
applied to density matices is generally called characteristic or 
generating function representation. Its properties as a repre- 
sentation for general operators have been extensively studied 
in |6| where due to its geometrical features it has been called 
the chord representation. We adhere to this nomenclature. In 



this paper we show the advantages and the simplicity of using 
this representation as a "noise" basis, where various models 
of noise are given a simple phase space interpretation. 

The paper is organized as follows. In Sec. Ill Al we briefly 
review some fundamental concepts about quantum open sys- 
tems and propose a two-stage scheme for the noisy propagator 
of density operators. In Sec. Ill Bl we describe noise models in 
phase space. Particularly, noise superoperators which are di- 
agonal in the chord or characteristic representation |6] or, as 
we show, whose Kraus operators are proportional to transla- 
tion operators in phase space. In Sec. [Ill] we find generaliza- 
tions of the depolarizing and phase damping channels |7, ^ 
for many qubits using the chord representation. Finally, in 
Sec.jiyjwe describe the diffusive noise recently proposed in 
1 9, 10] and show how it can be used to obtain the main part of 
the spectrum of the noisy propagator. 



II. A MODEL OF NOISY EVOLUTION 



A. Two-stage Superoperators 



Quantum systems coupled to a Markovian environment 
evolve according to the master equation Illlfl2ll 



dp 
It 
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LjPL] - -L]Ljp - -pL]L. 



(1) 
where Lj are the Lindblad operators. The first term on the 

right gives the unitary evolution. The second term represents 
interaction with an environment. Equation Q generates a so- 
lution in the form of a one parameter family (quantum dynam- 
ical semigroup 1 13]) of linear maps or superoperators \-t, such 
that the state at time t is given by 



pt = Lt(/5o). 



(2) 
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A general way of representing a superoperator is as a Kraus 
operator sum l\M . Given a superoperator L there exist a set 
of (Kraus) operators M^, such that 



^M^0Mt, 



(3) 



where, for L to be trace preserving, the identity ^ M^M^j, = 

/ must hold. Throughout the paper the symbol should be 
interpreted as 



MQM^p) = MpM^ 



(4) 



Not only does the Kraus form ensure that L(/3) is positive, 
for any density operator p, but it is also completely positive, 
meaning that tracing over any environment on which L acts 
trivially yields again a density operator [15]. 

In this paper instead of modeling environments and solving 
the master equation we focus on the properties of L for open 
quantum systems and thus propose models for noisy propa- 
gators. This shift of emphasis is analogous to the shift from 
hamiltonians to quantum maps in the study of unitary dynam- 
ics. 

We center our attention on discrete time systems 



Pn = L(/5„_i) = L"(/5o) 



(5) 



and propose a two-stage propagator consisting on the compo- 
sition of a unitary with a noisy evolution. After one unitary 
step a density matrix p becomes 



UpW 'i^' U(p) 



(6) 



The operator U is a unitary superoperator with trivial Kraus 
form. The interaction with a Markovian environment is mod- 
eled by a noise superoperator D^ depending on a parameter 
that quantifies its strength. The noisy one-step propagator is 
thus defined as 



L{p) = D. o [}{p). 



(7) 



Similar troes of two-stage schemes can be found in llgl lTolfm 
[T/, 18, 19]. There are several situations where this scheme 
appears naturally. One example is a kicked map in which the 
noise only acts between kicks. A billard inside a bath consti- 
tutes another example if we consider that the interaction with 
the walls is purely unitary while the free propagation is noisy. 
A quantum algorithm, supposed perfect, sent through a noisy 
channel is another example. In the case of a quantum algo- 
rithm that needs to be iterated, like Grover's search, the noisy 
part would be an effective interaction acting after each itera- 
tion. Finally a numerical solution of Lindblad's equation pro- 
ceeding in small time steps will naturally alternate between 
the unitary and the noisy propagation. 



B. Diagonal noise in the cliord representation 

Systems of K qubits are usually treated in a tensor product 
basis, the computational basis, defined as the tensor product 
of the eigenstates of (t^, { | 1) , |0) }. It is labeled by an integer 
\n) ranging from n = 0, ...iV — 1 where N = 2^ (throughout 
the paper, for many qubit systems, we call K the number of 
qubits and N the dimension of phase space). The coeficients 
of the binary expansion of n represents the state of each single 
qubit. As an operator basis, it is also customary to take the 



tensor product of Pauli operators, consisting of the N'^ = 4^ 
possible tensor products of 1, ax, o-y,o-z- 

An alternative approach, standard in the treatment of quan- 
tum maps 1 6 ] , and also recently utilized in the context of quan- 
tum information [4^ 20^ 211 is to treat ii'-qubit systems us- 
ing a phase space setting. In this approach the computational 
basis is assimilated to the position basis Bq = {\n),n = 
0, • • • , N — 1} (with periodic boundary conditions). There- 
fore the conjugate discrete momentum basis Bp = {\k),k = 
0, • • • , N — 1} is obtained by means of the discrete Fourier 
transform (DFT). The dimension N is then naturally related 
to an effective Planck constant by 



N = l/2Trh. 



(8) 



The natural operator basis in this context is constituted by 



the phase space translations 



l22r T(^q.p) defined as 



f(,,p)|fc) = cxp[-{27rt/N)q{k+p/2)]\k+p) (9) 
f(,,p)|n> = exp[{27n/N)pin + q/2)]\n + q), (10) 

with the group composition rule 



-'('3l,Pl)-'(92,P2) 



('3l+<?2,Pl+P2) *= I V-^-*^ 



where the phase in the exponential is the area of the triangle 
defined by the vertices (0,0), {qi,pi) and ((72,^2)- Transla- 
tions can be written in terms of compositions of translations 
in position and momentum. From Eq. Jilt . 






where 



C/9 = T, 
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yp aa J, 



(9:0), 



(12) 



(13) 



(o,p)- 



For the sake of simplicity, from now on we use a single (greek 
letter) index to represent the iV^ {<l,p) points of the discrete 
phase space, except where two indeces are needed explicitly. 
From Eq. M\\ it is easy to show that 



f^f = i. 



(14) 



These operators have matrix representation in TIn'^ or "Liou- 
ville" space, the space of N x N complex matrices with the 
Hilbert-Schmidt inner product 



(i,B)=Tr(its) 



(15) 



where A and B e T^Ara (we use "matrix" or "operator" to 
refer to elements in TYjv^ indistinctly, except where an explicit 
distinction is required). There are N"^ operators T(q.p) in Hn2 
and they form a complete orthogonal set since 



Tr(fif„,) = NS, 



(16) 



Therefore any operator A in 'Hn2 can be expanded in this 
unitary basis as 



^^^J2 "(")^" 



The c-number function (or symbol) 



1 



a(a) = -7^Tr(iTi) 



(17) 



(18) 



defines the chord representation [6] of A. The DFT of the 
translation operators (with the dimension extended to 2N, see 
1I23II ') are the generalized phase space point operators. Op- 
erators written in this basis constitue the Weyl or center 16] 
representation. The symbol of a density operator in the cen- 
ter representation is the well known discrete Wigner function, 
while in the chord representation it is also known as gener- 
ating or characteristic function. A description of the phase 
space point operators and the peculiar features of the discrete 
Wigner function can be found in f?|. 

A general superoperator D can be written in the basis of 
translations as 



Q,/3 



/3- 



(19) 



where the only requirement is that the matrix coeficients 
C(a,/3) be non-negative and C has unit trace. In this gen- 
eral setting the Kraus operators are linear superpositions of 
translations obtained by diagonalizing C{a, /3). In this paper 
we only consider superoperators where C{a, j3) ~ 5ajjC{a) 
is already in diagonal form in such a way that the Kraus oper- 
ators are simply proportional to the translations T^. Thus, the 
noise superoperator that we consider is explicitly defined by 
giving the Kraus operator sum form 



^a = ^Y.^M)T^QTl 



Trace preservation is achieved if 



Y,C.{a)lN^l. 



(20) 



(21) 



From Eq. \20\ we see that 0,^ is a convex sum of unitaries 
yielding a contracting superoperator Physically the action 
of Dct can be simply interpreted as performing an incoherent 
sum of all the possible translations T^ each with a probability 
Ca{ct)/N. In this way the choice of Ccr(a) determines differ- 
ent types of noise. The parameter <t is introduced to control 
the strength of the noise. 

One can think of this as a noise channel where the "errors" 
are unitary translations in phase space where they occur with 
probability Ca{a)/N . In this sense it constitutes a "very 
nice" 1 24] error basis with properties that are different from 
the more usual ones given by tensor products of Pauli matri- 
ces. 



One of the motivations for using noise in the form of 
Eq. ( 12m is that all the spectral properties are readily available. 
From Eq. (IIU it is clear that 



a 



N 

- C„{\)fx (22) 

where A is the wedge product XAfi = ^ip—vq, with A = (^, i^) 
and a = {q, p). Then the eigenfunctions of D^r are the transla- 
tion operators T\ and the corresponding eigenvalues are given 
by Ccr(A), the DFT of Ca{a). As Do- is diagonal in the chord 
representation ( I17t its action is quite simple. If p is expanded 
as 

then 



1 



(24) 



In the chord representation the action of Do- is simply to mod- 
ulate the elements px with Ca{X)- 

It is also of considerable interest to determine the spectral 
properties of the combined action of Do- with a unitary step U 
as in Eq. Q. In many instances (see sectionsllllland iryi a sig- 
nificant portion of the noise spectrum Cc (A) is zero or is con- 
tained within a small boundary of zero in the complex plane, 
with just a few isolated eigenvalues in the annular region be- 
tween it and the unit circle. This quasi-null subspace reduces 
the effective rank of the combined supeoperator Do o U and al- 
lows a very efficient calculation of the leading spectrum with 
diagonalizations of relatively small size (see Sec. lIV Al . 

III. GENERALIZED NOISE CHANNELS IN PHASE SPACE 

In this section we obtain many-qubit generalizations of the 
depolarizing and phase damping channels and write their su- 
peroperators in terms of translations in phase space. With 
a slightly different approach these generalizations have been 
studied in |25]. 

In Sec. IIIBI we proposed a noise superoperator which is 
diagonal in the chord representation and whose Kraus oper- 
ators are proportional to the translation operators. The main 
feature of this type of noise is the (diagonal) supermatrix of 
coeficients C{a) (or equivalently the spectrum which from 
Eq. (|23 is C{(3), i.e. the DFT of C{a)). We introduce a pa- 
rameter e which, if we think the noise is due to a coupling to 
an environment, it allows to continuously change from no cou- 
pling (e = 0) to full-strenght coupling (e = 1). Thus Eq. (I20> 
can be spht as a convex sum of two superoperators 

S, = (1 - e)to QT^ + ^Y. ^(")^" Tl (25) 

a=0 



Hence Se is also diagonal in the chord representation. Fol- 
lowing Eq. J22l i. the spectrum of S^ is given by 

E(/3) - (1 - e) + eC{(3) (26) 

(with 13 = (q,p) and q,p — Q^- ■ ■ ,N — 1). This means that 
the spectrum takes the constant value (1 — e) for all Tp plus 
an additional eC{P). Moreover C(0) = 1 for Sg to be unital. 

A. Depolarizing Channel 

The depolarizing channel for a single qubit, leaves it un- 
changed with probability (1 — e) and depolarizes it, which 
means that it leaves it in a completely mixed state, with proba- 
bility e. The Kraus operators for one-qubit depolarizing chan- 
nel Hi JJ are 



Mo 
Ml 

M2 

Ms 



VO^^l 



CTl 



0-2 



CT3 



2€ 

2(2^-1) 



<^1, 



2(2^-1)^2, 

2£ ;:;■ 

2(22-l)°^3- 



(27) 



If not stated explicitly matrices are written in the computa- 
tional basis {|0), |1)}. Notice that each A/^ is a constant mul- 
tiplied by the corresponding Pauli matrix ct^ (with &q = I), 
which are the generators of the SU{2) group and constitute 
an orthonormal basis of 7^22 . For K qubits a straightforward 
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Figure 1 : Wigner function representation of the action after one step 
of the depolarizing channel S , with N = 32 and e — 0.9. The 
left panel shows as initial state a cat state formed by the superpo- 
sition of two coherent states centered at (q,p) = (0.4, 0.25) and 
(0.6, 0.75). Long wavelength interference on the center corresponds 
to pure quantum interference while the othe short wavelength fringes 
correspond to interference with images generated by the periodicity 
of the torus. The uniform averaging over all phase space is clearly 
seen on the right panel. 



generalization of Eq. (|^) can be constructed with the Kraus 
operators given by 



where 7^ is the ^th element of the set of N^ — 1 generators of 
SU{N). They are hermitian operators which satisfy 



Ti-{7m} = 


(29a) 


Tr{%>} = 26 f,^. 


(29b) 


The normalized operators 




^"-> 


(30a) 


Q, - Ql- ^7. 


(30b) 



form a complete basis in Ti.i\j2 . We propose for the generalized 
depolarizing channel the following expresion 



,DC 



(1-6)70/+^ Y^Q^QQ, 

1-1=0 



def 



(l-e)/0/ + eS2 



(31) 



Using the fact that the set {Q^}^=£o spans the same subspace 
of traceles operators as {Ta/VN}a^o, it can be seen that S2 
in Eq. (TH is trace preserving. 

The computational basis, which we arbitrarily took to be the 
position basis in discrete phase space, defines a "canonical" 
orthonormal basis in Liouville space TYjv^ ■ The elements are 
the skew projectors (transition operators) 



with i, j = 0, 



p. 



P^j = N)(j|, (32) 

, iV — 1. Using Eq. dlOt it is easy to see that 

JV-l 



1 



f't p. . 



^ A y Tr 

q,p=0 
N-1 

= _ y g-*#p(«+j)y^ 



Ti 



ii^p) 



{i-i:P) ■ 



(33) 



p=0 



where l/v iV is added for normalization. For clarity, we used 
the two indices of T/^q.p) explicitly. 

Now, the generators 7^ can be written in the computational 
basis in terms of skew projectors as (see, for example, l26R ) 



7m 



{C/l2,C/l3,C/23, 



,^12, Vi3, V23, 

,Wn-i} 



with 



Ujk — Pjk + Pkj, 

Vjk = ijPjk-Pkj) 



Wi= - 



P22 + 



i(i+Tyv- 1 1 + ^22 - 
Pa - lPi+i,i+i) 



(34) 



(35) 



where l^j<A:<iV, 1^?<7V-1. Inserting Eqs. (|33} 
and ( I35> into Eq. J31> it can be transformed into 



Mq ^ y/ q^/ , 

JV£ -, 

2(Af2-l)7/ij 



M^ = 



(28) 



Sf = (l-6)f„0fot + ^ y f^QTl 

a=0 



(36) 



The physical interpretation is quite simple from a phase space 
point of view. With probability (1 — e) it leaves the state un- 
changed while with uniform probability e/N"^ it performs all 
the possible translations, thus averaging over all phase space 
with equal weight (except at the origin). In FIG.[Othe action 
of S on a superposition of two coherent states (represented 
by the Wigner function) is shown, with e = 0.9. The value of 
e is taken purposely large to make the effects of the averag- 
ing become apparent. Eventually, further averaging over the 
whole phase space leads to a completely depolarized state, i.e. 
the uniform state I /N . 

Equation \26\ gives the spectrum 



^DC 



w = (i-^) + ^e4^^^" 



/JV)/3Ao 



(1 - e) + eSp^o- 



(37) 



The spectrum takes the constant value (1 — e) for all T(q p) 
except for f'(o,o) where it is equal to 1 (FIGEI left). The (TV^ - 
l)-degeneracy makes the composition of this channel with any 
other superoperator trivial. The effect is to contract all the 
spectrum, except for a single 1, by a factor (1 — e). 

When e '^ 1 the spectrum lies within a small boundary of 
the origin leaving 1 as the unique non-zero eigenvalue. Phys- 
ically this means that after one step all th modes that are or- 
thogonal to the identity decay. Thus the only preferred state 
in this limit case is the uniform density I/N. 



B. Phase-Damping 

The phase-damping channel is widely used in the contexts 
of quantum to classical correspondence and quantum informa- 
tion because it provides a simple picture of how decoherence 
acts picking out preferred sets of states |7]. The Kraus opera- 
tors for the one-qubit phase damping channel are 



Mo ^^/{T^I, 
Ml =V^|0)(0|, 
Ma =V^|1)(1|. 



(38) 



It is instructive to see what happens after repeated action of 
the phase-damping channel. It is easy to check that after n 
steps, the initial density matrix po is 



(S 



PDCx 



'(/5o) 



Poo (1 - e)"/Ooi 

(1 - e)"pio Pii 



(39) 



So if e represents a decay rate, exponential decay of the non- 
diagonal terms occurs leaving the state in a completely mixed 
state. Therefore decoherence picks out the computational 
states as preferred basis. Unlike the depolarizing channel, 
phase damping only produces loss of coherence. 



We propose a generalized S 
The skew projectors Pij 



PDC 



for K qubits in terms of 



AT-l 



»PDC 



(l-e)/0/ + e 



2^ CijPij Q Pij- 

hj=0 



(40) 



To identify y^eCij Pij with Kraus operators then the coeffi- 
cients Cij must be real and positive. Moreover, in order to 
be trace preserving the identity ^ 
ing the same arguments of Sec 



Ci. 



1 must hold. Us- 
an expresion in terms 
of translation operators can be obtained. For example, if 

Ci-j = Si-j we get 



^PDC 



= (1- 



W-1 

e)^0^+4y"T'(o,p)0^/c 



N ^ 

p=0 



(0:P)- 



(41) 



The physical interpretation follows intuitively. With proba- 
bility (1 — e) it leaves the state unchanged, while it averages in 
a preferred direction (in this case the vertical lines q = const., 
p E [0, A^ — 1]). The channel picks out as pointer basis the 
position projectors Pa. This is not a surprise since we (arbi- 
trarily) took as computational basis, the position eigenstates 
\q). Had we chosen the momentum basis instead, then the 
averaging over lines p — const would lead to momentum as 
pointer basis. Following this train of thought we show that dif- 
ferent pointer bases can be selected by "diffusing" along any 
line in the phase space grid. We define a line on the N x N 
grid Gn of points {q,p) in phase space as 



nin2n3 == {(Q^P) ^ ^N] Tiip ^ ^12? + "-3, ^i, 7^2,^3 G N} 



L 

(42) 
The generalized phase damping channel on the hne imnana 
is given by 



(g,p)eL„j„2„3 

(l-e)/0/ + eD„,„,„3, (43) 



e 
R 



where Dn^n^ns is a superoperator like ( I20> . The integer R is 
the number of points contained in the line. If either rii or ria 
is odd then R=N. On the other hand, it is an integer power 
of 2 times N if both rii and 712 are even 0]. 

The eigenvalues corresponding to each translation T(^q,p) 
follow from the analysis in Sects. UrBl andlllll 

I]™C(g,p) = 1 - e(l - e-*(2V^)9n3/ni^^^^^^^^^) (44) 
for ni ^ 0, while if ni =0 and 712 7^ it is 

^^°^{q,p) = 1 - e(l - e'(2^/^)P"^/"M,^o). (45) 

This spectrum is displayed in FIG.|2](b) and (c). 

In FIG.|3]the effect of S„^„2n3 acting on a superposition of 
coherent states can be seen, for the case rii — l,n2 = —1 and 
n3 = 0. The right panel shows how the line (p = —q) that 
the channel picks, and averages over, appears. The only dif- 
ference with a line with non-zero 71,3 would be a vertical trans- 
lation of the line represented. Moreover, there is a progresive 
erasure of quantum interferences between the two coherent 
states. 

The action of this type of noise channel can be best under- 
stood if we delve deeper into the action of D„j„2„3. Using 




Figure 2: Complex plane plot of (a) The spectrum of the generalized depolarizing channel. Notice the A'''^ — 1 degeneracy on (1 — e, 0); (b) 
The spectrum for the Phase damping channel on the line with TV = 32, ni = 1, n2 = 2, m = 2. There are A'^^ — A'^ eigenvalues on (1 — e, 0) 
and TV (doubly degenerate) eigenvalues on the circle of radius e centered at (1 — e, 0); and (c) The spectrum for the Phase damping channel 
on the line with TV = 32, ni = 1, n2 = 0, na = 2. There are N'^ — N eigenvalues on (1 — e, 0) and TV eigenvalues equal to 1. 
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Figure 3: Wigner function representation of the action after one step 
of the phase damping on the line Si,-i,o, with TV = 64 and e = 0.85. 
The initial state is the same as in FIGQ The averaging over the line 
p = —q (mod 64) can be appreciated on the right panel. (Notice that 
the short wavelength interferences are due to the torus periodicity.) 



Eq. ( I12> it can be written as a composition of two superopera- 
tors 



D wria/ni 

nin20 O V , 

D I |-H3/n2 

On20 ° U ; 



ni =0 



(46) 



where the divisions in the exponents should be understood as 
the product with the inverse of an integer, modulo TV (when- 
ever such inverse exsists). Therefore V"^'"^ — y"3/"i 
(yt^ns/ni jg a unitary translation of na/rii in momentum and 
U-»3/«2 ^ ij--ns/n2 Q (-;7t)-n3/n2 is a Unitary translation of 
—ns/n2 in position. In other words, first apply a translation 
in position or momentum and then average on the line of slope 
n2/ni that contains the origin. 

This superoperator has interesting features when e ~ 1. 
Equations ( I44> and ( I45t show that there is an TV^ — R dimen- 
sional null subspace. The remaining non-zero eigenvalues are 
on the unit circle and, depending on the parity of rii , n2 and 
713, may exhibit double (or some integer power of 2) degen- 
eracy. The effect is that the averaging superoperator D„j„20 



truncates the unitray superoperator V"'^'"^ to a smaller sub- 
matrix which is still unitary (the same applies to Do„2 com- 



posed with U 



-ri3/n2 



When this noise channel acts on a unitary map (with e ~ 1) 
the full propagator is truncated to an i?-dimensional subspace 
spanned by translations on the chosen line. 



IV. DIFFUSIVE NOISE 

Although the effect of the previous two types of noise is 
clear in the context of many-qubit quantum algorithms, it is 
not intuitive to picture them as diffusion in phase space, be- 
cause they average over large distances (namely the whole 
space). On the other hand, although a purely diffusive noise 
is difficult to "see" acting on a system of qubits, it has a nat- 
ural interpretation in phase space. However, in spite of this 
differences, they can be treated as different types (depend- 
ing on a function C{a)) of the same class of noise given in 
Eq. ( I25l l. Following recent works j^ [13, flal we define as 
diffusive noise. 



D. = VaHTaor^ 



(47) 



which is a particular case of Eq. ( I25> with e = 1 and where 
Ca (a) is a periodic function narrowly peaked around a — 
and of approximate width a. For practical purposes we take 
Cct {a) to be a periodized Gaussian of half-width a and thus 
Cl{\) is also a Gaussian of half-width l/(27rfTTV) (FIG.|4| 
where A = (/i, ly)). This implements an incoherent sum of all 
the possible translations in phase space weighed by Ccr{a). 
To avoid a net drift C^ (a) should satisfy 



Ca{a) = Ccr{-a) 



(48) 



which makes 0^^ a Hermitian operator. 

In the context of quantum optics a similar type of noise for 
continuous phase space variables has been called Gaussian 
noise L27.1 . 
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Figure 4: Schematic plot of Ccr{fi,iy). The torus phase space is 
equivalent to a square with periodic boundary conditions. The points 
on the grid represent the states. The size of the dark region represents 
the effective rank of the truncated il(p',q'){p,q) (Eq. <50t ). Numeri- 
cally we set the white part equal to 0. 



A. Spectral properties of the noisy propagator 

In the examples of Sec.|III|the effect of noise acting on the 
unitary propagation was very simple. On the other hand, in 
the case of diffusive noise the numerical calculation is more 
involved and its nature depends strongly on the properties of 
the unitary map. In the context of quantum classical transition 
the spectrum of the noisy propagator, for chaotic systems, has 
gained interest because in the limits of ft ^ and no noise 
it can be related to purely classical properties like the Ruelle- 
PoUicott (RP) resonances |18, 28]. The RP resonances de- 
termine the long time decay of classical quantities 12911 and as 
was recently shown \W\ they determine the long time aproach 
to equilibrium of such quantities as the purity and fidelity for 
classicaly chaotic quantum systems with noise. Experimen- 
tally the appearance of RP resonances in quantum systems has 
been shown, for example, in I30ll . Clearly the part of the spec- 
trum that can be related to long time decays is the largest in 
modulus (and smaller than 1). 

The chord representation is specially suited for these cal- 
culations because since Ca (a) is a Gaussian of half- width a 
then the spectrum Cct(A) is also a Gaussian of complementary 
half- width l/(27rcr) which means that for some ranges of aN 
there is a large part of the spectrum of D^ very close to zero. 

The large (iV^ x iV^) supermatrix of the unitary step in the 
chord representation is 



^a' a — tl 



To^'UTM 



(49) 



Since the noise superoperator is diagonal in this representation 
its application is particularly simple and is given by 



-Ca'A — Crj{\ )ilA'A 



(50) 



Thus for A' 3> l/(criV) the matrix elements of Eq. \5Q\ are 
negligibly small. Therefore C^ (A) tends to kill the long chord 



components of the unitary step. A schematic representation 
of Ccr(A) is found in FIG. |3 The grid represents points in 
phase space. The whiter region gives the idea of which matrix 
elements can be neglected. The number of points inside the 
darker region represents the rank of the noisy superoperator 
reduced by the diffusion. 

In the numerical treatment we pick a threshold of values for 
Ccr(A). To do so we determine a range which is a coeficient 
a times the width of Ccr(A). Outside this range the quai-null 
spectrum is set exactly equal to zero. The resulting matrix is 
of a size of order 4(a/(27ro'A^))^ (see FIG. |3l. This proce- 
dure provides, up to a reasonable accuracy, the largest part of 
the spectrum of L. The calculation becomes more difficult as 
(T — > and to reobtain the unitary spectrum the full dimen- 
sion N"^ X A^^ is needed. For chaotic maps the spectrum is 
composed of the eigenvalue 1, corresponding to po = I /N, a 
small set of eigenvalues with modulus smaller than one, and 
many eigenvalues close to zero. The action of the noise is to 
set the latter ones exactly equal to zero. 

In FIG|5]we see computation of the spectrum for three dif- 
ferent sizes of truncation for a quantum map with diffusion. 
The leading part of the spectrum is unchanged while ther ap- 
pear slight in the less significant part. The stability of the 
spectrum up to orders of 10^^ can be clearly observed. 

Many different schemes have been used to calculate the 
leading spectrum of the noisy porpagator|9, 18, 31, 32. 13311 . 
The advantage over these methods is that although they use 
smaller matrices, the number of eigenvalues obtained is lim- 
ited by accuracy and diffucult error estimation. 

In 1 19] there is an analogous analysis for diffusive noise on 
a spherical phase space. 



V. CONCLUSIONS 

We provided a generalization of two well known noise 
channels in the context of phase space representations of 
quantum mechanics. Specifically we showed how the depo- 
larizing and phase damping channels have superoperator ex- 
pressions that are very simple to study using the chord repre- 
sentation. Moreover their spectral features, which can be cru- 
cial when composed with a map or an arbitrary algorithm are 
well determined. Some useful properties of these generalized 
noise channels composed with quantum maps and algorithms 
are currently being studied ]34]. 

The noises where formulated as a one parameter (e) family 
and their spectral properties where studied. In particular when 
the spectrum of the noise has a large null eigensubspace the 
convolution with a unitray map results in a truncated matrix 
whose effective rank depends on the rank of the noise. Thus 
reducing the computational requirements for systems with a 
large Hilbert space. This was shown to be useful when com- 
puting the leading spectrum of coarse grained propagators cor- 
responding to classsically chaotic quantum maps. The coarse 
graining was obtained by means of a diffusive noise and we 
showed that the leading spectrum is independent of the trun- 
cation. 
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Figure 5: Leading spectrum of L = Do-U, where 0^ is the difussive noise of Sec. IIVI (with a = 0.063) and U is the unitary superoperator 
for the quantum perturbed Arnold cat used in 13, 11011 (with A'^ = 100, k = 0.02). The dimension of the propagator after the truncation 



is dim = 4:{a/ (2n a)Y. (a) a = 2, dim = 4(2/27ro-)^ ^ 100; (b) a = 2.8 and dim = 4(2.8/27ro-)^ 
dim = 4(4.8/27r(T)^ f» 576. (Note that N^ x N^ = 10000.) If A,; is the i-th eigenvalue, then log \i = log(ri) 
the coordinates in the plots are (c^, — log(r)). 



f» 196; (c) a = 4.8 and 
i</f>i (where r-t — \Xi\) and 
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